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Abstract 

We investigate the cosmological perturbation dynamics for a universe consisting of pressureless 
>-' ' baryonic matter and a viscous fluid, the latter representing a unified model of the dark sector. In 

the homogeneous and isotropic background the total energy density of this mixture behaves as a 
generalized Chaplygin gas. The perturbations of this energy density are intrinsically nonadiabatic 
and source relative entropy perturbations. The resulting baryonic matter power spectrum is shown 
to be compatible with the 2dFGRS and SDSS (DR7) data. A joint statistical analysis, using also 
Hubblc-function and supernovae la data, shows that, different from other studies, there exists a 
> i maximum in the probability distribution for a negative present value , « -0.53 of the deceleration 

parameter. Moreover, while previous descriptions on the basis of generalized Chaplygin-gas models 

^o ■ . . . 

were incompatible with the matter power-spectrum data since they required a much too large amount 

of pressureless matter, the unified model presented here favors a matter content that is of the order 

f^*) | of the baryonic matter abundance suggested by big-bang nucleosynthesis. 



* E-mail: hipolito@ceunes.ufes.br 
' E-mail: velten@cce.ufes.br 

* E-mail: winfried.zimdahl@pq.cnpq.br 



I. INTRODUCTION 

Since 1998 a huge amount of data has been accumulated which directly or indirectly 
back up the conclusion, first obtained in []J, that our current Universe entered a phase 
of accelerated expansion. Direct support is provided by the luminosity-distance data of 
supernovae of type la (SNIa) [2] (but see also [3(), indirect support comes from the anisotropy 
spectrum of the cosmic microwave background radiation J] , from large-scale-structure data 
[5(], from the integrated Sachs- Wolfe effect [6], from baryonic acoustic oscillations [7(] and from 
gravitational lensing [8f]. Most current cosmological models rely on the assumption that the 
dynamics of the Universe is described by Einstein's general relativity and a material content 
that is dominated by two so far unknown components, pressureless dark matter (DM) and 
dark energy (DE), a substance equipped with a large negative pressure. For reviews of the 



actual situation see 9Hll| and references therein. The preferred model is the ACDM model 
which also plays the role of a reference model for alternative approaches to the DE problem. 
According to the interpretation of the data within this model, our Universe is dynamically 
dominated by a cosmological constant A which contributes more than 70% to the total cosmic 
energy budget. More than 20% are contributed by cold dark matter (CDM) and only about 
5% are in the form of conventional, baryonic matter. Because of the cosm olog ical constant 
problem in its different facets, including the coincidence problem (see, e.g., [12l.ll3||). a great 
deal of work was devoted to alternative approaches in which a similar dynamics as that of 
the ACDM model is reproduced with a time varying cosmological term, i.e., the cosmological 
constant is dynamized. Both DM and DE manifest themselves so far only through their 
gravitational interaction. This provides a motivation for approaches in which DM and DE 



appear as different manifestation of one sing! 



model and its different generalizations 
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e dark-sector component. The Chaplygin-gas 



23| realize this idea. Unified models of the dark 



sector of this type are attractive since one and the same component behaves as pressureless 
matter at high redshifts and as a cosmological constant in the long time limit. While the 
homogeneous and isotropic background dynamics for the (generalized) Chaplygin gas is well 



compatible and even slightly favored [24| by the data, the study of the perturbation dynamics 
resulted in problems which apparently ruled out all Chaplygin-gas type models except those 
that are observationally almost indistinguishable from the ACDM model [25J . The point here 
is that a generally finite adiabatic speed of sound in generalized Chaplygin gas (GCG) models 
predicts oscillations (or instabilities) in the power spectrum which are not observed. Also 



the analysis of the anisotropy spectrum of the cosmic microwave background disfavored these 
models [26l. |27||. except possibly for low values of the Hubble parameter [28(. To circumvent 
this problem, nonadiabatic perturbations were postulated and designed in a way to make the 
effective sound speed vanish [29. l3Qf| . But this amounts to an ad hoc procedure which leaves 
open the physical origin of nonadiabatic perturbations. There exists, however a different 



type of unified mode" 
was argued in 



31 



of the dark sector, namely viscous models of the cosmic medium. It 



321 ] . that a viscous pressure can play the role of an agent that drives the 



present acceleration of the Universe. The option of a viscosity-dominated late epoch of the 



Universe with accelerated expansion was already mentioned in [33j], long before the direct 
observational evidence through the SN la data. For more recent investigations see, e.g. [34 - 



391 ] and references therein. In the homogeneous and isotropic background viscous fluids share 



the same dynamics as GCGs 
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- 361] . But while perturbations in a (one-component) GCG 
are always adiabatic, viscous models of the dark sector are intrinsically nonadiabatic. In 
a recent paper we showed explicitly that, different from the Chaplygin-gas case, the power 
spectrum for viscous matter is well behaved and consistent with large-scale-structure data. In 
particular, it does not suffer from the mentioned oscillation problem 401 ] . On the other hand, 



what is observed in the redshift surveys is not the spectrum of the dark-matter distribution 
but the baryonic matter spectrum. Including a baryon component into the perturbation 
dynamics for a universe with a Chaplygin-gas dark sector, it turned out, that the mentioned 
oscillation within the dark component are not transferred to the baryons 
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23| . The 



baryonic matter power spectrum is well behaved and consistent with observation. Instead, 
there appears the new problem that the unified Chaplygin-gas scenario itself is disfavored 
by the data. It is only if the unified scenario with a fixed pressureless (supposedly) baryonic 
matter fraction of about 0.043 (according to the WMAP results) is imposed on the dynamics, 
that consistency with the data is obtained. If the pressureless matter fraction is left free, 
its best-fit value is much larger than the baryonic fraction. In fact it becomes even close 
to unity, leaving only a small percentage for the Chaplygin gas, thus invalidating the entire 
scenario. In other words, a Chaplygin-gas-based unified model of the dark sector is difficult 
to reconcile with observations. One may ask now, whether the status of unified models can 
again be remedied by replacing the Chaplygin gas by a viscous fluid. It is exactly this question 
that we are going to investigate in the present paper. It is our purpose to study cosmological 
perturbations for a two-component model of baryons and a viscous fluid, where the latter 
represents a one-component description of the dark sector. We shall show that such type 



of unified model is not only consistent for a fixed fraction of the baryons but also for the 
case that the matter fraction is left free. Our analysis demonstrates that the statistically 
preferred value for the abundance of pressureless matter is compatible with the mentioned 
baryon fraction 0.043 that follows from the synthesis of light elements. 

The contents of the paper is as follows: in section[II]we establish our two-component model 
of a viscous dark component and baryons and discuss its background dynamics. Section 
IIIII is devoted to the perturbation dynamics of this mixture. Subsection IIIIAI considers 
the nonadiabatic total energy-density perturbations, subsection IIIIBI presents a dynamical 
equation for the relative entropy perturbations and in subsection llll Cl we obtain the fractional 
baryonic energy-density perturbations which are shown to be adiabatic at high redshifts. A 
numerical integration and tests against data from the matter power spectrum, the Hubble 
function H(z) and SNIa are given in section HV] which also contains a statistical analysis of 
the validity of the viscous unified model itself. Finally, section [V] summarizes and discusses 
our main results. 

II. THE TWO-COMPONENT MODEL 

The cosmic medium is assumed to be describable by an energy-momentum tensor 

T ik = pu l u k + ph ik , h ik = g ik + u l u k , (1) 

which splits into a matter part T$ and viscous fluid part Ty, 

T ik = T ik + T ik ; ( 2 ) 

with 

Tl) = p M ui lU k M + p M (g ik + t4-t&) > T v = Pvu v u v + Pv (g ik + u v u k v ^ , (3) 

where the subscript "M" stands for matter and the subscript "V" stands for viscous. The 
total cosmic fluid is characterized by a four velocity u m while u l M represents the four velocity 
of the matter part and u v represents the four velocity of the viscous fluid. Energy-momentum 
conservation is supposed to hold separately for each of the components, 

T ik . i =Ty k . i = Q => T%=0. (4) 

In particular, the energy balances are 

PM,iU Z M + U L Mi (p M + Pm) = , PV,iUv + u V:i (PV + Pv) = (5) 



and 

p ti u i + rf i (p + p) = , (6) 

where (up to first order) p = pu + Pv and p = pu + Pv- In general, the four velocities of the 
components are different. We shall assume, however, that they coincide in the homogeneous 
and isotropic zeroth order, 

u M = Uy = u l (background) . (7) 

Difference will be important only at the perturbative level. 

Let the matter be pressureless and the viscous fluid be characterized by a bulk viscous pressure 

Pv, 

p M = , py = p = -CO , (8) 

where £ = const and © = vH: i is the fluid expansion. Under this condition the total pressure 
coincides with the pressure of the viscous component. The total background energy density 
is p = pm + Pv, where 

p v + 3H (p v + pv) = , pm + 3-ff pm = =>- pm = PMoa~ 3 . (9) 

The total energy balance is p+3H (p + p) = 0. In the homogeneous and isotropic background 
one has O = 3H, where H is the Hubble rate. If, moreover, the background is spatially flat, 
the Friedmann equation 3H 2 = 8irGp implies oc p 1 ' 2 , such that p = —Q(2AttG) ' p 1 ' 2 . 
This coincides with the special case a = — \ for the equation of state p = — -^ of a gener- 



alized Chaplygin gas, if we identify A = Q y/2An G. In terms of the present value qo of the 
deceleration parameter q = — 1 — -^ the total energy density can be written as [40(] 
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l-2q + 2(l + q )a- 3 2] , => — = - [l - 2g + 2 (1 + g ) tHl , (10) 



2 H 1 

where po and Hq denote the present values of p and H, respectively. Since pM = Pmo&~ 3 , we 
have pv = p — /Qa/oQ 3 - These relations show that it is the total energy density that behaves 
as a GCG, not the component V. This type of unified model differs from unified models 
in which the total energy density is the sum of a GCG and a baryon component. Only if 
the baryon component is ignored, both descriptions coincide. For the total equation of state 
parameter we obtain 

P 1 - 2g 



P l-2g + 2(l + g )o" 



(11) 



Consequently, in the homogeneous and isotropic background, a generalized Chaplygin gas 
with a = —1/2 can be seen as a unified description of the cosmic medium, consisting of a 
separately conserved matter component and a bulk viscous fluid with C, = const, where the 
latter itself represents a unified model of the dark sector. 

III. PERTURBATIONS 

A. Nonadiabatic perturbations of the total density 

The system is characterized by the equations of state ([8]). It is expedient to emphasize 
that we have neither an equation of state py = pv(pv) nor an equation of state p = p(p). It is 
only in the spatially flat background when, via Friedmann's equation, the relation p = — £0 
reduces to p oc —p 1 ' 2 and the corresponding energy density coincides with the energy density 
of a GCG. Neither the component V nor the system as a whole are adiabatic. Because of 
p = — C©j the pressure perturbation is p = — £0, where a hat on top of the symbol denotes the 
(first-order) perturbation of the corresponding quantity. The nonadiabaticity of the system 
as a whole is characterized by 

* P ' -P- P ,D = 3H^( P --t), (12) 



p+p pp+p p p \p Q 

where we have introduced the abbreviations 

P = -?— , D = -?— . (13) 

p+p P + P 

The quantity (|12|) is governed by the dynamics of the total energy-density perturbation p and 
by the perturbations of the expansion scalar, which is also a quantity that characterizes the 
system as a whole. The behavior of these quantities is described by the energy-momentum 
conservation for the entire system and by the Raychaudhuri equation, respectively. Both of 
these equations are coupled to each other. The remarkable point is that these quantities and, 
consequently, the total energy density perturbation, are independent of the two-component 
structure of the medium. The reason is the direct relation p = — £0 between the pressure 
perturbations and the perturbations of the expansion scalar. This is different from pertur- 
bations in a two-component system where each of the components is adiabatic on its own. It 
will turn out that the total energy-density perturbations are characterized by a homogeneous 
second-order differential equation. These perturbations, which are intrinsically nonadiabatic, 
then act as source terms in the evolution equation for the relative entropy perturbations. The 

6 



perturbations in the baryon component are obtained as a combination of the total and the 
relative entropy perturbations. 

The general line element for scalar perturbations is 

ds 2 = - (1 + 20) dt 2 + 2a 2 F a dtdx a + a 2 [(1 - 2if>) 5 a[S + 2E >a$ ] dx^dx^ . (14) 

Since g mn u m u n = — 1 and also g mn u 1 2u r ^ = — 1, it follows that 

no = u = u M = Uy = —<j> and a u^ + a F jfi = u^ = v^ . (15) 

The last relation defines the quantity v which will be used to introduce gauge invariant 
quantities on comoving (v = 0) hyper surf aces. Similarly, one defines the corresponding 
quantities vm and vy for the components. These different velocity potentials are related by 

. Pv +Pv , \ a Pu , , n „s 

vm = v H [vm — W) an d vy = v [vm — vy) . (16) 

p+p P+P 

We also introduce the quantity 

x = a 2 (E-F) . (17) 



The combination v + x 1S g au g e invariant. It is convenient to describe the perturbation 
dynamics in terms of gauge invariant quantities which represent perturbations on comoving 
hypersurfaces, indicated by a superscript c. These are defined as 

P c _ P . 6 C _ 6 p c _p 

— = - + v, - 7 - = ? - + v, — = - + v. (18) 

p p G 6 p p 

For the fractional quantities we introduce the abbreviations 

D c = J— , P c = -?— . (19) 

p+p P+P 

In our case we have 

(20) 
In terms of the comoving quantities the total energy and momentum balances may be com- 



P _ e ff_ e^ 

P ~ e P ~ e 



bined into (cf. [40j]) 



D c -3H^D c + e c = . (21) 

P 



The expansion scalar Q is governed by the Raychaudhuri equation, 

9 + ^6 2 + 2 (a 2 - u?) - ii% + 4vr G (p + 3p) = . (22) 



Up to first order the perturbed Raychaudhuri equation can be written in the form 

. (23) 



<§> c + 2HO c + \AP C + ^H 2 D c 



It is through the Raychaudhuri equation that the pressure gradient comes into play: 

P 



7/? e 



1 v 2 
P c = — ^D c 
2 1P 2 



3-fpH' 



-D c 



(24) 



where 7 = 1 + | . The pressure perturbation consists of a term which is proportional to the 
total energy-density perturbations D c (notice that the factor in front of D c is positive), but 
additionally of a term proportional to the time derivative D c of D c . The relation between 
pressure perturbations P c and energy perturbations D c is no longer simply algebraic, equiv- 
alent to a (given) sound-speed parameter as a factor relating the two. The relation between 
them becomes part of the dynamics. In a sense, P c is no longer a "local" function of D c but 



it is a function of the derivative D c as well 4l(. This is equivalent to p = p(p,p). It is only 



for the background pressure that the familiar dependence p = p{p) is retained. As already 
mentioned, the two-component structure of the medium is not relevant here. 
Introducing now 



5 = ~/D c 



L 
p 



(25) 



and changing from the variable t to a, Eqs. (|2ip and (|23p may be combined to yield the 
second-order equation 



5" + f(a)S' + g(a)d = 0, 



(26) 



where 5' = -j- and the coefficients / and g are 



/(«) 



and 



9(a) 



3 p 1 p k 2 
2 ~% ~ 3tp H 2 a 2 



3 15 v 9 p 2 1 p 2 k 2 



(27) 



(28) 



respective' 



case in [4C ] . 



2 2/o 2 p 2 -yp 2 H 2 a 2 _ 
y. Equation (|26p coincides with the corresponding equation for the one-component 



B. Relative entropy perturbations 

Alternatively to relation (|12p , the deviation from adiabaticity in a two-component system 
with components M and V is 



P 



P P 



P + P P P + P 



pc 



p 



D c 



PV +PV 



PV 



PV PV 



P + P \Pv+Pv PvPv+Pv 



+ 



Pm (fiv + Pv) PV 
{p + pf Pv 

Solving this for the nonadiabatic part of component V yields 

P + P 



Pv 



PV +PV 



PM 
PM 



PV 



PV PV 



PV + PV PV PV +PV PV + PV 



P c - -D c 



3H-— I PM Pv 



P P \PM PV 



The perturbed energy balances for the components (A = M, V) are 



PA 



+ 3H 



PA 



PA PA 



3ip + -T (Ava + A X ) = 



(29) 



(30) 



(31) 



^PA+PAj \PA+PA PA PA + PA 

Obviously, the combination (|30p enters the energy balance of the viscous component. Sub 
tracting the balance of fluid M from the balance of fluid V and using (|30p it follows that 



PV 


_ Pm\ 

Pm) 


+ 3H< 


P + P 
^Pv +PV 


P c - 


P 


_ 3H PPM_ [PM 
P P \PM 

+—A(v v - 


_Pv\ 
Pv) _ 

- vm) = 


PV +PV 



. (32) 

To deal with the term that contains the difference vy — vm of the velocity potentials of the 
components, we implement the momentum balances which imply (A = M, V) 

PA , PA 



+ 



-VA + VA + 



. 



PA +PA PA+ PA 
With pm = 0, the definition for P c in (|19|) and with f|16j) we arrive at 



p + p p p M , v 

P -3H- {vm-vv) 



{v v - v M )' = 

PV +PV PPV + PV 

Introducing relative entropy perturbations by the usual definition 

_ PM pv 

omv = 



(33) 



(34) 



(35) 



PM PV + PV 

differentiating equation (j32[) and combining the result with equation (]34p and with (j32[) again, 
we obtain the inhomogeneous second-order equation 



S'vm + r ( a ) s 'vM + s(a)S VM = c(a)5 f + d(a)8 



(36) 



with the coefficients 



r(a) = — 
a 



s(a) 



[3_3P_ 3 


P PM 


2 2p ppv+p 


3 p PM 
a 2 p pv +p 


_1 + 7- 



(37) 



(38) 



and 



c[a) 



d(a) 
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P 



P 



.1 Pv +P 



l + f +! + -*■ 



A- 2 



9tf 2 a 2 



27a 2 py + p 



l- p - 



1+-2- 

2p 



2^(l + ^ 



A: 2 



(39) 



(40) 



p J \ 'ipj p \ jpj 9H 2 a 2 
The set of equations (j36[) and (|26[) contains the entire perturbation dynamics of the system. 
At first, the homogeneous Eq. fl26j) for 5 has to be solved. Subsequently, once 5 is known, 
Eq. (I36p determines the relative entropy perturbations. 



C. Baryon density perturbations 



The quantity relevant for the observations is the fractional perturbation 5m — „ of the 
energy density of the baryons. This quantity is obtained from the total fractional density 5, 
determined through (l26j) . and the relative entropy perturbations Svm, determined through 
(ESD, by 



5m = - 

7 



r Pv +P Q 

bvM 



with 



p v + p _ 2 (1 + q ) a" 3 / 2 [1 - 2q + 2 (1 + q ) a" 3 / 2 ] - 9ft M0 a- 3 



[l-2go + 2(l + go)a- 3/2 ]' 



(41) 



(42) 



ttG 



where we have introduced the present value of the matter fraction Qmo = %jjiPM0- Assuming 



3H,i 



Hq to be given, the free parameters of the system are qo and Qm ■ 

At early times, i.e. for small scale factors a <C 1, the equation 
form 

s " + l 6 '-i? s=0 ' ia « r > 

independent of go an d for all scales. The solutions of (|43p are 

5 (a < 1) = cia + c 2 a~ 3/2 , 



has the asymptotic 



(43) 



(44) 
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where c\ and ci are integration constants. The nonadiabatic contributions to the total density 



perturbations are negligible at high redshifts 40]. 



For n«l the coefficients s(a), c(a) and d(a) in ([36]) become negligible and r(a) — > |. 
Eq. ([36]) then reduces to 

S VM + ^S' VM = 0, (a«l). (45) 

It has the solution Svm = const = 0. From the definition (|35p we find that at high redshifts 

S MV = ™--PL, (««!), (46) 

PM PV 

since -Jf- <C 1 under this condition. Consequently, there are neither nonadiabatic contribu- 
tions to the total energy-density fluctuations nor relative entropy perturbations and we have 
purely adiabatic perturbations 5 m = i5 at a < 1. This allows us to relate our model to the 
ACDM model at early times. We shall use the fact that the matter power spectrum for the 
ACDM model is well fitted by the BBKS transfer function [4J]. Integrating the ACDM model 
back from today to a distant past, say z = 1.000, we obtain the shape of the transfer function 
at that moment. The spectrum determined in this way is then used as initial condition for 
our viscous model. This procedure is similar to that described in more detail in references 



43, 



4J. 



IV. STATISTICAL ANALYSIS 



To estimate the free parameters of our model we perform a Bayesian analysis and construct 
the corresponding probability distribution functions. At first we consider large-scale-structure 
data from the 2dFGRS [45[| and SDSS DR7 [46] programs. The matter power spectrum is 



defined by 



Pk = \5 M ,k\ 2 > (47) 



where 5m,h is the Fourier component of the density contrast 6m- Generally, for a set of 
free parameters {p}, the agreement between the theoretical prediction and observations is 
assessed by minimizing the quantity 

*,._ 1 [Pf(p)-^( P) ] 2 

x (p)-jv7 2^ ~2 ' ( 48 ) 

i % 

where Nf means the number of degrees of freedom in the analysis. The quantities P\ and 
pobs are fa e theoretical and the observed values, respectively, of the power spectrum and &{ 

11 



denotes the error for the data point i. With the help of x 2 we then construct the probability 
density function (PDF) 

p = S e -4 Ei , (49) 

where B is a normalization constant. 

To test our model against the observed power-spectra data we consider the following 
two situations, (i) We assume the matter component to be entirely baryonic with a fraction 
^Afo = 0.043 as suggested by the WMAP data. Fixing also Hq = 72, a value favored by these 
data as well, the only remaining free parameter is go- This will provide us with information 
about the preferred value(s) of go for the unified dark-sector model, (ii) We leave the matter 
fraction free, thus admitting that the matter component is not only made up by the baryons. 
This is equivalent to allow for a separate DM component in addition to the contribution 
effectively accounted for by the viscous fluid. This additional freedom is used to test our 
unified model of the dark sector itself. The unified model can be regarded as favored by 
the data if the PDF for the matter fraction is large around the value that characterizes the 
baryon fraction. If, on the other hand, the PDF is largest at a substantially higher value, 
the unified model has to be regarded as disfavored. The PDF for case (i) is shown in figure 
[TJ We obtain two regions with high probability for go, one of them with a pronounced peak 
around qo ~ —0.53, implying accelerated expansion. The other one, which is of the same 
hight, has qo > and is compatible with an Einstein-de Sitter universe. The appearance of 
a maximum of the PDF in the region go < is neither observed in Chaplygin-gas scenarios 
nor in our previously studied one-component viscous model [40j]. The difference to the latter 
might appear surprising since go characterizes the system as a whole and the addition of a 
small fraction of baryons should, at the first glance, not have a large impact on the total 
dynamics. However, it is not the background dynamics that counts here. In the present 
case the PDF for go is inferred from power-spectra data that are related to the fluctuations 
5m, while in the one-component model these data were related to the fluctuations 5 of the 
total energy density. As relation (|4ip shows, 5m and 5 may be very different in general. The 
appearance of a maximum for go < means, that the results of a first-order analysis may 
well be compatible with the results for the background, which imply go < 0. We consider 
this an advantage over Chaplygin-gas models, for which there has always been a tension 
between the results in the background and those on the perturbative level 



22, 



23J. Figure[2] 



(Figure [3]) shows the theoretically obtained spectrum for various negative (positive) values of 
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qo together with the power-spectrum data points. To better illustrate the relation between 
the predictions of the model and the observations, two different normalization wave numbers, 
k n = 0.034/iMpc~ 1 and k n = 0.185/iMpc -1 , have been chosen, but our statistical results do 
not depend on a specific normalization. 

In order to break the degeneracy between the high-probability regions in Fig. Q] we include 

n 

information from the Constitution set of SNIa (see the last reference in 21) and from the 

I 

recent H(z) data from [47]. The results from the joint analysis with SNIa data are shown in 
Fig. [J while Fig. [5] depicts the corresponding PDF for qo, based on the joint analysis with 
the H(z) data. In Fig. 0]the total \ 2 is calculated from x 2 = x\dFGRS + Xsdss + XsNia' m 
Fig. \5\ from x 2 = X 2 2dFG RS + X 2 SD SS + Xh- 

For case (ii) we have both qo and £Imo as f ree parameters. The results of the statistical 
analysis are shown in figure [6l Most importantly, there is a high probability for small values 
of the matter fraction Ojv/O) including the WMAP value f^A/o = 0.043. According to our 
previously mentioned criteria this means, the unified viscous model is indeed preferred by 
the data. This is in striking contrast to unified Chaplygin-gas models which have high 



probabilities close to Qmo = 1, t 
of dark matter and dark energy 



ius apparently invalidating the idea of a unified description 



22 



23|. 
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FIG. 1: One-dimensional PDF for qo resulting from the 2dFGRS data (solid curve) and from the SDSS DR7 
data (dashed curve). The right picture is an amplification of the peak in the region go < 0. 
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FIG. 2: Power spectra (PS) normalized at k n = 0.034hMpc 1 (left panels) and at 0.185hMpc 1 (right panels) 
for different negative values of qo- The top panels compare the PS with the 2dFGRS data, the bottom panels 
with the SDSS DR7 data. 



V. DISCUSSION AND CONCLUSIONS 

It is well known that, in a homogeneous, isotropic and spatially flat background universe, 
there exists an equivalence between viscous and generalized Chaplygin-gas models for a 
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PS with k n =0.034 hMpcT 1 (2dFGRS) 
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FIG. 3: Power spectra (PS) normalized at k n = 0.034/iMpc 1 (left panels) and at 0.185/iMpc 1 (right panels) 
for different positive values of go. The top panels compare the PS with the 2dFGRS data, the bottom panels 
with the SDSS DR7 data. 
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FIG. 4: Left panel: Hubble diagram, center panel: PDF for go, based on a joint analysis of PS and SN data. 
The right panel magnifies the maximum for go < 0. 



15 



2dFGRS + SDSS DR7 + H(z) data 
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FIG. 5: Left panel: Hubble parameter as a function of the redshift for different values of go. Right panel: 
one-dimensional PDF for qo, based on a joint analysis of PS and H(z) data. 
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FIG. 6: PDF for the pressureless component Amo (left) and for the deceleration parameter qo (center), using 
the 2dFGRS data (solid lines) and the SDSS DR7 data (dashed lines). The right picture is a normalized 
amplification of the peak for qo < in the central panel. 
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361 ] . The cosmic substratum at the 



unified description of the cosmological dark sector 
present time is then approximated as a mixture of one of these dark components and baryons. 
The novel approach presented here is based on the fact that also the two-component system 
of a bulk viscous fluid and a separately conserved baryon component behaves in the back- 
ground as a generalized Chaplygin gas with a = —\. The total energy-density perturbations, 
however, are intrinsically nonadiabatic and coincide with those of a one-component viscous 
fluid, investigated in earlier work J40(]. While the baryon component may be considered 
dynamically negligible in the background, the situation is different on the perturbative level, 
since the observed matter agglomerations are related to baryonic density fluctuations. These 
fluctuations are obtained from a combination of the said nonadiabatic total energy density 
perturbations and relative entropy perturbations in the two-component system where the 
former source the latter. The observed matter-power spectrum is well reproduced. There do 
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not appear oscillations or instabilities which have plagued adiabatic Chaplygin-gas models 

jrove the findings of a previous one-component analysis in which 



251 ] . Our present results im 



no baryons were included 40j]. At first, the probability distribution for the deceleration 
parameter has a maximum at qo ~ —0.53 which partially removes the degeneracy of previous 
studies which, taken at face value, were incompatible with an accelerated expansion and 
thus in obvious tension with results for the background. Perhaps still more important is 
the test of the unified model itself. Many investigations on approaches with a unified dark 
sector fix the pressureless matter component to be that of the favored (by the WMAP data) 
baryon fraction and then check whether or not the resulting dynamics can reproduce the 
observations. This corresponds to the strategy (i) in the previous section. But this is only 
part of the story since it does not say anything on how probable the division of the total 
cosmic substratum into roughly 96% of a dark substance and roughly 4% of pressureless 
matter is. To decide this question, one has to consider the pressureless matter fraction as a 
free parameter and to find out which abundance is actually favored by the data. Our analysis 
in point (ii) was devoted to this task and revealed that the matter fraction probability is 
indeed highest for values smaller than roughly 8%. This is a result in favor of the unified 
viscous model. We recall that a corresponding analysis for a Chaplygin gas results in values 
close to unity 



231 ] which seems to rule out such type of approaches. The present viscous 
model, on the other hand, remains an option for a unified description of the dark sector, at 
least as far as the matter power spectrum is concerned. 



Acknowledgement: We thank FAPES and CNPq (Brazil) for financial support. WSHR is 
grateful to FAPESB and to S. Carneiro for interesting discussions and for their hospitality 
during his stay at the Institute of Physics at UFBA. 



[1] A.G. Riess et al, Astron. J. 116, 1009 (1998); S. Perlmutter et al, Astrophys. J. 517, 565 (1999). 
[2] J.L. Tonry et al, Astrophys. J. 594, 1 (2003); M.V. John, Astrophys. J. 614, 1 (2004); P. Astier 



et al, J. Astron. Astrophys. 447, 31 (2006); A.G. Riess et al, |astro-ph/0611572[ D. Rubin et 
aL larXiv:0807. 11081 M. Hicken et al, Astrophys.J. 700, 1097 f2009). [arX"iv:0901 .48041 

[3] S. Sarkar, Gen. Relativ. Gravit. 40, 269 (2008). 

[4] S. Hanany et al, Astrophys. J. Lett. 545, L5 2000; C.B. Netterfield et al, Astrophys. J. 571, 



17 



604 (2002); E. Komatsu et al., Astrophys. J. Suppl. 180, 330 f2009h larX"iv:0803.0547l 
[5] M. Colless et al, Mon. Not. R. Astron. Soc. 328, 1039 (2001); M. Tegmark et al, Phys. Rev. D 
69, 103501 (2004); S. Cole et al., Mon. Not. R. Astron. Soc. 362, 505 (2005); V. Springel, C.S. 
Frenk, and S.M.D. White, Nature (London) 440, 1137 (2006). 
[6] S. Boughn and R. Chrittenden, Nature (London) 427, 45 (2004); P. Vielva, E. Martinez- 
Gonzalez, and M. Tucci, Mon. Not. R. Astron. Soc. 365, 891 (2006). 



[7] D.J. Eisenstein et al, Ap.J. 633, 560 (2005), |arXiv:astro-ph/0501171 



[8] C.R. Contaldi, H. Hoekstra, and A. Lewis, Phys. Rev. Lett. 90, 221303 (2003). 

[9] E.J. Copeland, M. Sami and S. Tsujikawa, Int.J.Mod.Phys.D15, 1753 (2006). 

[10] T. Padmanabhan, Gen. Relativ. Gravit. 40, 529 (2008). 

[11] R. Durrer and R. Maartens, Gen. Relativ. Gravit. 40, 301 (2008). 



[12] N. Straumann, astro-ph/0203330 



[13] T. Padmanabhan, Phys.Rept. 380, 235 (2003); [heplh /0212290 

[14] A. Yu. Kamenshchik, U. Moschella and V. Pasquier, Phys.Lett. B511, 265 (2001). 

[15] J.C. Fabris, S.V.B. Gongalves e P.E. de Souza, Gen. Rel. Grav. 34, 53 (2002). 

[16] M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev. D66, 043507 (2002). 

[17] L.M.G. Bcca, P.P. Avelino, J.P.M. dc Carvalho and C.J.A.P. Martins, Phys. Rev. 67, 101301(R) 

(2003). 
[18] N. Bilic, G. P. Tupper, and R. D. Viollier, Phys. Lett. B535, 17 (2002); N. Bilic, R. J. Lindebaum, 

G. P. Tupper, and R. D. Viollier, J. Cosmol. Astropart. Phys. 11 (2004) 008; M. C. Bcnto, O. 

Bertolami, and A. A. Sen, Phys. Rev. D 70, 083519 (2004); R. R. R. Reis, M. Makler, and I. 



Waga, Phys. Rev. D 69, 101301 (2004); R. A. Sussman, |arXiyi080L3324; N. Bilic, G. P. Tupper, 

and R. D. Viollier, arXiv: 0809.0375. 
[19] D. Carturan and F. Finelli, Phys. Rev. D 68, 103501 (2003). 
[20] W. Zimdahl and J.C. Fabris, Class. Quant. Grav. 22, 4311 (2005). 
[21] V. Gorini, A.Y. Kamenshchik, U. Moschella, O. F. Piatella, and A. A. Starobinsky, J. Cosmol. 

Astropart. Phys. 02. 016 (2008). 
[22] J.C. Fabris, S.V.B. Gongalves, H.E.S. Velten and W. Zimdahl, Phys. Rev. D 78, 103523 (2008). 
[23] J.C. Fabris, H.E.S. Velten and W. Zimdahl, Phys. Rev. D 81, 087303 (2010). 
[24] R. Colistete Jr, J. C. Fabris, S.V.B. Gongalves and P.E. de Souza, Int. J. Mod. Phys. D13, 669 

(2004); R. Colistete Jr., J. C. Fabris and S.V.B. Gongalves, Int. J. Mod. Phys. D14, 775 (2005); R. 

Colistete Jr. and J. C. Fabris, Class. Quant. Grav. 22, 2813 (2005); R. Colistete Jr. and R. Giostri, 



BETOCS using the 157 gold SNe la Data : Hubble is not humble, arXiv:astro-ph/0610916 



[25] H.B. Sandvik, M. Tegmark, M. Zaldariaga and I. Waga, Phys. Rev. D 69, 123524 (2004). 

[26] L. Amendola, F. Finelli, C. Burigana and D. Carturan, JCAP 0307, 005 (2003). 

[27] R. Bean and O. Dore, Phys. Rev. D 68, 023515 (2003). 

[28] M.C. Bento, O. Bertolami and A.A. Sen, Phys. Lett. B 575, 172 (2003). 



18 



[29] R.R.R. Reis, I. Waga, M.O. Calvao e S.E. Joras, Phys. Rev. D68, 061302 (2003). 

[30] L. Amcndola, I. Waga and F. Finelli, JCAP 0511, 009 (2005). 

[31] W. Zimdahl, D.J. Schwarz, A.B. Balakin, and D. Pavon, Phys. Rev. D 64, 063501 (2001). 

[32] A.B. Balakin, D. Pavon, D.J. Schwarz, and W. Zimdahl, NJP 5, 85.1 (2003). 

[33] T. Padmanabhan and S. M. Chitre, Phys. Lett. A 120, 433 (1987). 

[34] J.C. Fabris, S.V.B. Goncalves and R. de Sa Ribeiro, Gen. Rel. Grav. 38, 495 (2006). 

[35] M. Szydlowski and O. Hrycyna, Ann.Phys. 322, 2745 (2007). 

[36] R. Colistete Jr., J.C. Fabris, J. Tossa and W. Zimdahl, Phys. Rev. D76, 103516 (2007). 

[37] A. Avelino and U. Nucamendi, JCAP 0904 (2009) 006. 

[38] B. Li and J.D. Barrow, Phys. Rev. D79, 103521 (2009). 

[39] A. Avelino and U. Nucamendi, JCAP 1008 (2010) 006. 

[40] W.S. Hipolito-Ricaldi, H.E.S. Veltcn and W. Zimdahl, JCAP 06 (2009) 016. 

[41] W. Zimdahl, Int. J. Mod. Phys. D (IJMPD) 17, 651 (2008) (|arXiv:0705.213ip . 

[42] J.M. Bardeen, J.R. Bond, N. Kaiser and A.S. Szalay, Astrophys. J. 304, 15 (1986); J. Martin, 

A. Riazuelo and M. Sakellariadou, Phys. Rev. D61, 083518 (2000). 

[43] J.C. Fabris, I.L. Shapiro and J. Sola, JCAP 0702, (2007) 016. 

[44] H.A. Borges, S. Carneiro, J.C. Fabris and C. Pigozzo, Phys. Rev. D77, 043513 (2008). 

[45] S. Cole et al., Mon. Not. R. Astron. Soc. 362, 505 (2005). 

[46] B.A. Reid et al., Mon. Not. R. Astron. Soc. 404, 60 (2010). 

[47] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski and S. A. Stanford, JCAP 1002, 008 (2010). 



19 



